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Abstract 



Ph ' The Iwasawa decomposition of the symplectic group is used to write down 

closed- form expressions for the noncompact part of the group. The results 
make possible the computation of Lyapunov exponents using the method de- 
. scribed in |l| in systems with more than two degrees of freedom. 

o 

In these Proceedings [|1| Habib and Ryne propose a method for computing Lyapunov 
exponents in cases where the linearized dynamics is Hamiltonian. The method requires a 



closed-form expression for the noncompact part of the symplectic group, denoted in [[I]] by 
Sp(2n), in the defining representation. (Other common notations for the symplectic group 



O 

are Sp(n,~R) and Sp(2n,M,).) The authors provide such an expression for Sp{2) and Sp(4) 

O 

and state that they are working on generalizations to Sp(2n). In this note I point out that 
suitable closed-form expressions for the noncompact part of Sp(2n), n > 1, can be found 
rather easily using the Iwasawa decomposition of this group pfl. 

Denoting the set of matrices of the form JS C and JS a (Eq. (10) of U) by K and P, I 
first demonstrate that they represent the Cartan decomposition of the symplectic Lie algebra 
sp(2n). In the calculations that follow use is made of the Killing form for sp(2n) which is 
given by (Xi,X 2 ) = (2n + 2)Tt(XiX 2 ) for any X\ % G sp(2n) 0]. Then, since any element of 
K can be written in the form JS+SJ where S is a symmetric 2n x 2n matrix and J is denned 
in JI]], and any element of P in the form JS — SJ 0], it can be verified that [K, K] C K, 
[K, P] C P, [P, P] C K, (X, X) < for X ^ 0, X E K, and (Y, Y) > for Y ^ 0, Y G P. 
Thus sp(2n) = K + P is indeed the Cartan decomposition of sp(2n). K is the maximal 

1 



compact subalgebra, which in this case is the unitary Lie algebra u(n). To compute the 
closed-form expression for the noncompact factor of Sp(2n) displayed in Eq. (11) of I 
begin by writing down the Iwasawa decomposition of sp(2n), sp(2n) = u(n) + Hp + N. Hp 
is an abelian and N a nilpotent Lie algebra for which the generators can be taken as 



1 < i < n for Hp; 



(la) 



and 



■m 



ij &i,n+j &j,n+i) 



1 < i,j < n 



(lb) 



1 < hj < n i i < j for N. 



(lc) 



Here eij is the usual Weyl basis for 2n x 2n matrices, (e^) rs = 5i r 5j S . (The method for 
computing the Iwasawa decomposition for any semisimple real Lie algebra is given in ||.) 

An arbitrary element of Sp(2n) can now be written as e x e z e w , where X G u(n), Z G 
Hp, W G N ||. Since H P is abelian, the exponential series in e z is easily evaluated, whereas 
the nilpotency of N converts the infinite series in e w into a polynomial of order 2n — 1. The 
factor e z e w is thus 

2n-l 



e z e w = diag(Ai,...,A n ,A 1 1 ,...,A n 1 )(^/+ — ( ^ a^m^ + fey ray)' 



(2) 



s=l i,j i<j 

,-1 



The parameters ay and fey are in R and Ai,...A n > 0. The quantity of interest in [|TJ, 
MM, takes the form 



MM = [e z e w )e z e w . 



(3) 



As an example of the use of Eq. (|2]) I write down the matrix display for Sp(4): 

( 1 &12 2an - |a 2 2^ 2 a 12 + 022&12 \ 

1 a i2 - a 2 2&i2 2a 2 2 
1 

V0 -612 1 



diag(Ai, A 2 , A/, A 2 X ) 



(4) 
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